Elliptic flow of thermal dileptons in relativistic nuclear collisions 
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We calculate the transverse momentum and invariant mass dependence of elliptic flow of thermal 
dileptons for Au+Au collisions at the Relativistic Heavy Ion Collider. The system is described 
using hydrodynamics, with the assumption of formation of a thermalized quark-gluon plasma at 
some early time, followed by cooling through expansion, hadronization and decoupling. Dileptons 
are emitted throughout the expansion history: by annihilation of quarks and anti-quarks in the 
early quark-gluon plasma stage and through a set of hadronic reactions during the late hadronic 
stage. The resulting differential elliptic flow exhibits a rich structure, with different dilepton mass 
windows providing access to different stages of the expansion history. Elliptic flow measurements 
for dileptons, combined with those of hadrons and direct photons, are a powerful tool for mapping 
the time-evolution of heavy-ion collisions. 

PACS numbers: 25.75.-q,12.38.Mh 



I. INTRODUCTION 

Evidence has been mounting for some time that col- 
lisions between heavy nuclei at ultra-relativistic energies 
lead to formation of quark-gluon plasma (QGP). Strong 
confirmations arise from recent observations at the Rel- 
ativistic Heavy Ion Collider (RHIC) of large anisotropic 
flow of all hadronic species P, [3, H, 0] and of a suppres- 
sion of high-pT hadrons due to parton energy loss in the 
dense medium 0, Q . Signatures of direct photon emis- 
sion [7, 8, 9, as well as preliminary results on (excess) 
dilepton production in such collisions, indicative of a hot 
early state, have also started emerging 

The radial and elliptic flow of hadrons observed in rel- 
ativistic heavy-ion collisions at RHIC strongly suggest 
that the quark-gluon plasma (QGP) created in these col- 
lisions acts like a strongly coupled plasma with almost 
perfect liquid behaviour [l2| . This conclusion is based 
on the successful prediction of the hadron momentum 
distributions, in particular of their anisotropics in non- 
central collisions, by dynamical calculations which treat 
the expanding QGP as an ideal fluid 0, i, [H, [H, [H, HI . 
While no other equally successful model exists, one has 
to remain conscious of the fact that the new "perfect liq- 
uidity" paradigm is based on a model back-extrapolation 
of the measured data to the early stages of the colli- 
sion which are not directly accessible with hadronic ob- 
servables. There are strong arguments that this back- 
extrapolation is fairly unique [12| and hence that the 
above-mentioned qualitative conclusion is robust. On a 
quantitative level, however, the extraction from exper- 
imental data of the (small) QGP viscosity is presently 
hampered not only by the unavailability of consistent hy- 
drodynamic codes for viscous relativistic fluids, but even 
more by uncertainties about the hydrodynamic effects of 
changes in the equation of state of the QGP matter 
and about details of the initial conditions at the begin- 
ning of the hydrodynamic expansion stage [H, [l^ . 



The next generation of RHIC experiments will aim 
at obtaining more quantitative information about the 
properties of the QGP and the subsequent hot hadronic 
matter, as well as on the process of hadronization itself. 
While more precise hadronic flow measurements, in par- 
ticular of elliptic flow in non-central collisions, will fur- 
ther constrain the hydrodynamic evolution models and 
the QGP equation of state, it will be invaluable to have 
additionally measurements on electromagnetic radiation, 
i.e. direct photons and dileptons which are penetrating 
probes [l9j, and access more directly the very early ex- 
pansion stages 20]. Elliptic flow is generated very early, 
via the transformation of the initial spatial eccentricity of 
the nuclear overlap region into momentum anisotropies 
through the action of azimuthally anisotropic pressure 
gradients. With the passage of time, the pressure gra- 
dients equalize, and the growth of elliptic flow shuts it- 
self off [2l|. Photons and dileptons, which escape from 
the expanding fireball without re-interaction, will be 
able to probe specifically this early stage where the flow 
anisotropy first develops. Measurements of the elliptic 
flow of thermal photons and lepton pairs could thus pro- 
vide particularly clean constraints on the QGP equation 
of state, with fewer ambiguities arising from theoretical 
uncertainties about the hadronization process and late- 
stage hadronic rescattering dynamics than in the case of 
hadronic flow measurements [16]. 

The azimuthal asymmetry of photons from jet-plasma 
interactions has been examined recently (2^ . and so has 
the elliptic flow of thermal photons [23|, calculated with 
a hydrodynamic model for the expansion of the colli- 
sion fireball. Contrary to the monotonic rise with pT of 
the hydrodynamically predicted differential elliptic flow 
V2{pt) of hadrons, the thermal photon elliptic flow was 
shown to peak at transverse momenta of 1-2 GeV/c, de- 
creasing at higher px as a reflection of the decreasing 
hydrodynamic flow anisotropies as one goes backwards 
in the expansion history of the collision fireball. In ad- 
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dition, the photon eUiptic flow was seen to exhibit an 
interesting peak-valley structure at low transverse mo- 
menta, pt ^ 0.4 — 0.5GeV/c, reflecting the dominance 
of different types of hadronic photon production mecha- 
nisms below and above this pt value 1 23] . In the present 
paper we complement the work of [23] with a study of 
thermal dilepton (i.e. virtual photon) elliptic flow, where 
the invariant mass of the lepton pair (photon virtuality) 
provides an additional continuously tunable parameter. 
As we will show, the combined dependences of the ellip- 
tic flow of dileptons on their transverse momentum pT 
and invariant mass M provide a very rich landscape of 
structures which can be used to set observational win- 
dows on specific stages of the fireball expansion. Future 
studies will further complement the analysis presented 
here by directly measuring the space-time structure of 
the early collision fireball with two-photon correlations. 
By combining the momentum structure of the (virtual) 
photon emitting source from flow measurements with its 
space-time structure from Hanbury Brown- Twiss correla- 
tions one should be able to provide complete experimen- 
tal characterization of the phase-space structure of the 
QGP during the early stages of the flreball expansion. 



II. DILEPTON EMISSION RATES AND 
SPECTRA 



Dileptons, like photons are emitted from every sta ge o f 
a heavy- ion collision, from the pre-equilibrium stage [24| . 
the quark-gluon fluid [25, 26, 27, 28, 29 30, 31, S^, and 
the late hadronic matter [331, [3^, |35j, ISfii |37|, |38| . An in- 
teresting outcome of dilepton emission studies has been 
the conflrmation of medium modifications of the spectral 
properties of p mesons [13, [s^- Correlated charm and 
bottom decay provide another important source of dilep- 
tons, which will help estimate energy loss and the elliptic 
flow of heavy quarks fiol] . It may be possible in principle 
to isolate them by determination of the decay vertices. 

In the present paper we focus on thermal emission of 
dileptons from the QGP and hadronic phases in a col- 
lectively expanding flreball. The dilepton spectrum is 
obtained [34| by integrating the thermal emission rate 
over the space-time history of the system. As a result, 
high-M dileptons arise mostly from the hot early stage 
where hydrodynamic flow is weak but the spatial eccen- 
tricity of the source is large, whereas dileptons of lower 
mass are emitted when the temperature is low, the flow is 
strong and anisotropic, but the spatial eccentricity of the 
flreball has mostly disappeared. These generic expecta- 
tions are, however, modulated by the fact that dilepton 
emission from the late hadronic phase is characterized by 
strong vector meson resonance peaks in the dilepton mass 
spectrum which leads to strong variations of the relative 
weight of dilepton emission from the QGP and hadron gas 
phases on and off these resonances. This leads to very 
strong and interesting structures in the pr-integrated el- 
liptic flow V2{M) as a function of dilepton mass which 



can be exploited for differential fireball microscopy. 
The dilepton momentum spectrum can be written as 
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where the quantity inside the square brackets indicates 
the thermal emission rates from the QGP or hadronic 
matter, with the thermal Boltzmann factor extracted for 
emphasis. We parametrize the dilepton (virtual photon) 
4-momentum as 

p^ = {Mt cosh y, Pt cos 0, pr sin 0, Mt sinh Y) (2) 

and the flow velocity of the fireball fluid as 

= ^T(^cos\vri,Vx(x,y,T),Vy{x,y,T),smhri^ (3) 

with 
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(We assume boost-invariant longitudinal expansion ^47].) 
We use coordinates t, x, y, rj with volume element d*x = 
Tdrdxdydrj where r = (t^— z^)^/^ is the longitudinal 
proper time and r]= tanh~^(z/i) is the space-time ra- 
pidity. The dilepton momentum is parametrized by its 
rapidity Y, its transverse momentum px = {pi + Py)^^'^, 
and its azimuthal emission angle (j). 

The dilepton energy in the local fluid rest frame, which 
enters the Boltzmann factor in the thermal emission rate 
in the combination p-u/T, is thus given by 

p ■ u ^ [Mt cosh(y— 77) — ptVt cos{(j)—(t>v)] , (5) 

where Mt = (Af^ + PtY^^ is the transverse mass 
and M the invariant mass of the dilepton, and 
(pv = tan^^^Vy/vx) is the azimuthal angle of the trans- 
verse flow vector. This expression shows that the az- 
imuthal anisotropy (0-dependence) of the dilepton spec- 
trum, conventionally characterized by its Fourier coeffi- 
cients Vn (where only even n contribute at Y = 0), 
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l + 2v2{M,pT,b) cos(20) 



(6) 



is controlled by an interplay between the collective fiow 
anisotropy and the geometric deformation of the tem- 
perature field T{x,y,T) (see [2^ for more details). It 
vanishes in the absence of radial fiow, vt = 0. 

We use the boost invariant hydrodynamic code AZHY- 
DRO [13] which has been used extensively and suc- 
cessfully to describe midrapidity hadron production at 
RHIC and more recently to predict the elliptic fiow of 
thermal photons [1^. We study Au-|-Au collisions at 
-v/s = 200 A GeV and use the same initial conditions as in 
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[23 |. namely an initial time ro = 0.2fm/c for the begin- 
ning of the hydrodynamic stage, with an initial peak en- 
tropy density in central collisions of sq — 351 fm~'^, corre- 
sponding to a peak initial temperature of Tq — 520 MeV. 
Although the matter may not yet be completely thermal- 
ized at such an early time, we choose tq so small in or- 
der to also account at least partially for pre-equilibrium 
dilepton production at very early times ,24, 44j]. This 
contribution to the dilepton spectrum is important at 
large M and and large px , and it will suppress there 
because very little transverse flow develops before about 
0.5fm/c. The initial transverse entropy density profile is 
computed from a Glauber model assuming 75% wounded 
nucleon and 25% binary collision scaling of the initial 
entropy production [l3|. This reproduces the collision 
centrality dependence of hadron production in Au-|-Au 
collisions at RHIC [ll|4|. 

We assume that a thermally and chemically equili- 
brated plasma is formed at the initial time, and use 
the Born term for the production of dileptons from the 
QGP. It is known [sol, [sH that for thermal dileptons 
the corrections to the Born term are ~ (27ras/3)T^/Af ^ 
and thus decrease rapidly with increasing dilepton mass. 
For the production of dileptons from the hadronic re- 
actions, we consider the comprehensive set of hadronic 
reactions analyzed by Kvasnikova et al. [38| . which are 
known to correctly reproduce the spectral density mea- 
sured in e+e~ scattering. In [sl] these rates are con- 
veniently parametrized in the form Fcff{M) exp{~E/T). 
They have been shown to provide a good description of 
intermediate mass dilepton data from Pb+Pb collisions 
at the SPS ^3^. However, since these rates are used here 
for lower invariant masses, some comments are in order. 
For instance, it is known that vector meson spectral den- 
sities are modified in the medium |19|. and such effects 
are not included in the rates of Ref. [33| ■ Still, as long as 
they do not produce severe distortions of the transverse 
momentum spectra, the dilepton elliptic flow signal eval- 
uated here should be relatively robust. 

For the equation of state we use EOSQ 
which matches a free quark-gluon gas (QGP) to a 
chemically equilibrated hadron resonance gas (HG) 
by a Maxwell construction at critical temperature 
Tc = 164MeV, with energy densities Cq = 1.6 GeV/fm^ 
and Ch = 0-45 GeV/fm'^ in the QGP and HG subphases 
at this temperature. Hadron freeze-out is assumed to 
happen at ef = 0.075 GeV/fm^ [H. 

Figure [T] shows the production of thermal dileptons 
from the quark matter and the hadronic matter for cen- 
tral collisions. Similar figures emerge at other impact pa- 
rameters and help us understand the relative importance 
of hadronic and quark-matter contributions at a given M. 
The QGP contribution (dotted line) is seen to completely 
dominate the mass spectrum of thermal dilepton radia- 
tion for dilepton masses above 1 GeV, except for the 
meson peak. Hadronic radiation (dashed line) dominates 
for 0.5 GeV < M < 1 GeV. The red dash-dotted lines 
show post- freeze-out decay contributions from p, u, and 
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FIG. 1: (Color online) The mass spectrum of thermal dilep- 
tons from a hydrodynamical simulation of central 200 A GeV 
Au+Au collisions (& = 0) . The quark and hadronic matter 
contributions are shown separately. See text for details. 

(j) mesons; they are not included in the solid line for the 
total hydrodynamic spectrum. The sharply peaked solid 
black line under the lu peak (which is included in the to- 
tal hydrodynamic spectrum) shows the contribution from 
thermally distributed lu mesons in the hadronic phase 
before freeze-out. Since the hadronic dilepton emission 
rates given in [sl] do not include 3-body scattering chan- 
nels proceeding through the uj meson, we include the uj 
explicitly as a thermally distributed particle species, with 
its vacuum spectral function, integrating its density mul- 
tiplied with the standard partial width for dilepton decay 
over the space-time volume of the hadronic phase con- 
tained inside the hadronic freeze-out surface. We see that 
post-freeze-out decays of lu and (j) produce many more 
dileptons at the respective vector meson peaks than col- 
lisions and decays before hadronic freeze-out; this reflects 
the relatively long lifetimes of these two vector mesons 
(23.4 and 44.5 fm/c, respectively). 

We pause here for a moment to note that the rela- 
tive contributions of the emission from the hadronic and 
quark matter stages control the overall V2 for the dilep- 
tons. Thus we shall see in the following that the elliptic 
flow of dileptons with masses near those of the p, w, or 
(j) mesons is decided by the V2 of the radiation from the 
hadronic matter. The elliptic flow parameter for dilep- 
tons with M^l GeV, on the other hand, is dominated 
by quark matter radiation (see Fig. [5]). To measure this 
experimentally would obviously be very valuable. 

In Figures [2] and [3] we show the transverse momentum 
spectra and py-dependent elliptic flow V2{pt) for dilep- 
tons with invariant masses M — nip (Fig. [J) and M — 
(Fig. [3]), from semiperipheral Au-|-Au collisions at im- 
pact parameter 6 = 7 fm. The bottom panels of these 
figures show that the elliptic flow of hadronic dileptons, 
W2(HM), is large and basically agrees with the elliptic flow 
of hadrons of the same mass {v2{p) and V2{4>), respec- 
tively), which are emitted from the hadronic freeze-out 
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FIG. 2: (Color online) The transverse momentum spec- 
trum (upper panel) and differential elliptic flow V2{pt) (lower 
panel) of thermal dileptons with invariant mass M — mp . 
Dashed and dotted lines show quark matter (QGP,QM) and 
hadron matter (Had,HM) contributions, respectively. In the 
bottom panel we also show for comparison the elliptic flow of 
p mesons emitted from the hadronic decoupling surface. 



surface. The somewhat smaller differential elliptic flow 
V2{pt) for the hadrons (compared to that of dileptons 
with the same mass) arises from their somewhat flatter 
Pt spectrum (not shown). Since the hadrons decouple 
later than the average hadronic dilepton, their spectra 
are boosted by somewhat larger radial flow, and flatter 
spectra result in smaller differential elliptic flow V2{pt) 
(see last paper in Ref. Q). The elliptic flow of thermal 
dileptons from the QGP phase, W2(QM), is much smaller 
and shows the same decrease at large pr that we already 
observed for thermal photons [1^ . It also decreases with 
increasing dilepton mass. These features reflect the small 
flow anisotropics during the early QGP stage, especially 
for dileptons with large pT and M whose production re- 
quires high temperatures and is therefore peaked at very 
early times. One should keep in mind, however, that 
the present analysis does not include lepton pairs pro- 
duced via the Drell-Yan and jet fragmentation processes, 
nor those produced through jet-plasma interactions 43]. 
These can contribute their own characteristic azimuthal 




1 2 3 4 5 

Pt (GeV) 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1. 

■ Thermal Dileptons 


1 1 1 1 1 


: Au+Au@ 200 AGeV 




- b=7 fm 




rM=m^ . •• 




VaCHM).;' ^^^^ - 

















1 2 3 4 5 

Pt (GeV) 

FIG. 3: (Color online) Same as Fig. O but for dileptons with 
invariant M = m^, and with the elliptic flow of (/> mesons 
from the hadronic decoupling surface shown in the bottom 
panel for comparison. 



asymmetries [221. ^ consistent inclusion of these non- 
thermal contributions is part of a study in progress. 

The total dilepton flow V2 is seen to follow that of the 
hadronic dileptons, W2(HM), up to rather large transverse 
momenta of about 2 GeV/c, before peaking near 3 GeV/c 
and then decreasing at even larger px- This is due to 
the strong dominance of hadronic dilepton emission near 
the vector meson peaks in the dilepton mass spectrum of 
Fig. [1] As seen in the top panels of Figs. [2] and [3l the 
dilepton pT spectrum is dominated by hadronic emission 
over almost the entire shown pT range, with QGP ra- 
diation winning out only for transverse momenta above 
~ 4GeV/c. Note that a possible medium-induced sup- 
pression of the hadronic dilepton contribution by ~ 50% 
near the p mass would decrease the net V2 there by only 
10%, because it is a weighted average. Still, the absence 
of non-thermal sources such as those enumerated at the 
end of the preceding paragraph must be kept in mind 
and should be corrected before attempting a quantita- 
tive comparison with future dilepton flow data. 

In Figure S] we show that the success of the hydrody- 
namic model in describing all measured hadron px spec- 
tra, at least up to pt — 2GeV/c, carries over to that of 
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FIG. 4: Transverse momentum spectrum of mesons, com- 
paring the theoretical prediction from the hydrodynamic 
model with PHENIX data Ref. El. 



the (j) meson, recently measured by the PHENIX Collab- 
oration [46| via reconstruction of mesons from K~^K~^ 
pairs. It can thus be hoped that the predicted spec- 
trum and elliptic flow of dileptons with the same mass, 
as shown in Fig. [31 will be similarly reliable. 
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FIG. 5: (Color online) Same as Fig. [1 for M = 2 GeV. 

For dileptons with invariant mass M = 2 GeV the pt 
spectrum is completely dominated by emission from the 
QGP (Fig. O top panel). Correspondingly, the total 



dilepton elliptic flow V2 closely follows that of the quark 
matter contribution, W2(QM), giving us a direct measure- 
ment of the flow anisotropy developed during the quark 
phase. Its decrease at high px again reflects the decreas- 
ing flow anisotropics at earlier and earlier times. The 
total elliptic flow is small, peaking near pT — 2 GeV/c at 
a value of about 0.5%, even though the elliptic flow of the 
(tiny) hadronic contribution of the same mass is about 
20 times larger. 

The upper panel of Figure [5] shows the p^-integrated 
elliptic flow as a function of dilepton mass, V2{M). The 
solid line gives the total elliptic flow of all dileptons (ex- 
cluding, however, post-freeze-out decay dileptons), while 
the dashed and dotted lines show the elliptic flow of the 
QGP and hadronic dileptons separately. Note that the 
latter show similar qualitative behaviour as functions of 
M and of pt- while the elliptic flow of hadronic dilep- 
tons increases monotonically with M and pt, the quark 
matter dileptons exhibit elliptic flow that first rises, then 
peaks and finally decreases with increasing M and/or 
Pt [47l|. For comparison, that panel also indicates the 
PT-integrated elliptic flow values for a variety of stable 
hadron species and hadronic resonances emitted from the 
hadronic decoupling surface (see also the last paper in 
Ref. [1]). For large invariant masses the pT-integrated el- 
liptic flow of hadronic dileptons is seen to approach that 
of hadrons with the same mass. This reflects the fact 
that the homogeneity regions [131 for massive particles 
are small since their thermal wavelength decreases like 
1 / \/M, and therefore all particles of that mass (hadrons 
or virtual photons) feel approximately the same flow. 
Low-mass hadronic dileptons, on the other hand, are 
emitted from much larger homogeneity regions and thus 
have an appreciable chance of being emitted significantly 
earlier than the corresponding hadrons of the same mass; 
on average they thus carry less radial and elliptic flow. 

We note that the additional buildup of radial and (to a 
lesser extent) elliptic flow during the hadronic stage has 
opposite effects on the pT-integrated and pr-differential 
elliptic flows: While the pr-integrated elliptic flow in- 
creases, the additional broadening of the single-particle 
spectrum from extra radial flow shifts the weight of the 
flow anisotropy to larger py, so at fixed px the elliptic flow 
decreases (at least in the low-pT region where V2{pt) is 
a rising function oi pr)- This follows directly from the 
discussion presented in the last paper of Ref. , and it 
is clearly seen in the figures given in the middle and the 
lower panels of Fig. [HI The middle panel shows the (arbi- 
trarily normalized) pt distributions of dileptons having 
masses equal to uip and m^, respectively, compared to 
those of the p and mesons themselves. The hadron pT- 
spectra are clearly seen to decrease more slowly with pT^ 
as stated above. (We stress that the continued hydrody- 
namical growth of radial flow during the hadronic matter 
stage is essential for this pattern; if one neglects trans- 
verse flow one arrives at the opposite conclusion [Flj.) 
The lower panel in Fig. [H] then demonstrates that this 
flow-induced flattening of the hadronic p^-spectra leads 
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to a systematic lowering of the difFerential elliptic flow 
V2{pt) for these hadrons when compared to dileptons 
of similar masses. We shall see later in Fig. [7] (lower 
panel) that this is accompanied by larger values for the 
PT'-integrated flow paramters for p and (f) mesons com- 
pared to that of dileptons having similar masses. 

As seen in the upper panel of Fig. [51 the elliptic flow 
of the total dilepton spectrum shows dramatic structure 
as a function of dilepton mass, oscillating between the 
quark and hadronic matter limits. This reflects the rel- 
ative weight of the QM and HM contributions to the 



dilepton mass spectrum: Close to the vector meson res- 
onance peaks, dilepton emission is strongly dominated 
by emission from the hadron matter phase, and the solid 
line almost reaches the dotted curve for V2(HM). At low 
dilepton masses, between the p and (f) mesons, and for 
dilepton masses 1.5 GeV, quark matter radiation dom- 
inates the dilepton spectrum, resulting in small elliptic 
flow coefficients (of order 1% or less) for the total dilep- 
ton spectrum, consistent with the pure quark matter con- 
tribution W2(QM) (dashed hue). By setting appropriate 
invariant mass windows on the dileptons one can thus 
peek into the early QGP phase and study the beginning 
of elliptic flow buildup. 
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FIG. 7: (Color online) Top panel: Dilepton mass dependence 
of the pT-integrated elliptic flow parameter for thermal dilep- 
tons, for a variety of impact parameters (fe = 0, 2, 4, 6, 8, 10 
fm from bottom to top (solid curves) and 6 = 1, 3, 5, 7, 9 
fm (dashed curves)). Note that V2 for 6 = lOfm is marginally 
smaller than for 6=9 fm. Bottom panel: Impact parameter 
dependence of the pT-integrated elliptic flow parameter scaled 
by the initial spatial eccentricity e, for a variety of hadrons 
and dileptons with different masses. 

Figure [7] shows the impact parameter dependence of 
the pT-integrated elliptic flow parameter for thermal 
dileptons, as a function of dilepton mass M. The top 
panel shows that, as the impact parameter (and thus the 

spatial eccentricity e — of the initial nuclear over- 

lap region) increases, the elliptic flow V2 increases, too. 
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as naively expected. In the bottom panel we therefore 
study the eccentricity-scaled elliptic flow, V2/e, as a func- 
tion of collision centrality, for several hadron species and 
dileptons with a variety of different masses. Even though 
the dileptons are, on average, emitted earlier (i.e. with 
less radial flow) than the hadrons, we see no dramatic 
differences between the impact parameter dependences 
of hadronic and dilepton V2- Both are almost indepen- 
dent of collision centrality, except for the small drop of 
V2/e for dileptons with masses M < 1 GeV at large im- 
pact parameters 6^ 8fm. This reaffirms the general 
understanding [53 | that for ideal fluids the ratio V2 /e 
reflects the effective stiffness (sound speed) of the fireball 
medium's equation of state, averaged over the expansion 
history. The changing weight between QGP and hadronic 
dilepton emission for different dilepton masses is seen to 
only affect the magnitude of their elliptic flow, but not 
its collision centrality (in-) dependence. [Note that the uj 
contribution was not included in Fig. [T] The also omit- 
ted Drell-Yan contribution is expected to emerge only at 
larger invariant masses than those shown in Fig. [7| [53|.] 

III. CONCLUSIONS 

We have presented a first hydrodynamic calculation 
of elliptic flow of thermal dileptons emitted from ultra- 
relativistic heavy-ion collisions at RHIC energies. The 
azimuthal flow parameter exhibits a rich structure as a 
function of transverse momentum and invariant mass. 
When combined with photon elliptic flow measurements 
[23| . this yields a new versatile and potentially very 
powerful probe of the fireball dynamics at RHIC and 
LHC, complementary to the already well-studied flow 
anisotropics in the hadronic sector. vUpt) and V2^{M) 
exhibit rich structures which reflect the interplay of dif- 
ferent emission processes, opening a window on detailed 
and differential information from a variety of different 



stages of the fireball expansion. The elliptic flow of pho- 
tons emitted from the late hadronic stage was seen [2^ 
to track the V2 of the emitting hadrons, suggesting the 
possibility of subtracting the hadronic emission contribu- 
tions from the total photon signal in order to isolate and 
study in greater detail the elliptic flow of early QGP pho- 
tons. Similarly, we saw in the present study that at low 
{M < 0.5 GeV) and high (Af > 1.5 GeV) invariant masses 
dilepton emission is completely dominated by emission 
from the early quark matter, without need for subtract- 
ing hadronic contributions. However, at low M experi- 
mental backgrounds are large while at high M thermal 
dilepton yields are small. Clean elliptic flow measure- 
ments of the early QGP stage are therefore hard. 

Obviously, when compared to hadrons, electromag- 
netic probes suffer from their relatively small produc- 
tion cross sections, so these measurements require dedi- 
cation. However, the first glimpse offered here suggests 
that photon and dilepton elliptic flow have the potential 
of turning into powerful and highly discriminating tools 
for heavy-ion phenomenology and for differential studies 
of the fireball expansion. Extensive further theoretical 
investigations are therefore warranted, in particular with 
respect to the inclusion of viscous effects at intermediate 
and high transverse momenta. 
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